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In this work, a planar phased array of circular waveguides arranged in 
an equilateral triangular grid is analyzed. The boundary value problem is 
first formulated rather generally in terms of a vector two dimensional integral 
equation for an array of elements that are arranged in a doubly periodic grid 
along two skewed (nonorthogonal) coordinates. Dielectric plugs, covers, and 
loading, as well as thin irises at the aperture, are accounted for in the 
formulation. Numerical solutions are obtained by using the Ritz-Galerkin 
method to solve the integral equation. Excellent agreement with experimental 
measurements using a waveguide simulator is observed. The existence of 
forced surface wave phenomena in equilateral triangular grid arrays and 
their strong dependence upon the mode of excitation is also demonstrated. 
These phenomena are shown to exist at isolated points in the scan co- 
ordinates. Reflection characteristics as well as the polarization characteristics 
of the radiation pattern are illustrated at selected planes of scan for both 
linear and circular polarization excitation. 

I. INTRODUCTION 

The requirements of modern radar and communication systems have 
stimulated considerable activity in the design and use of phased 
array antennas. To date, the design information required for their 
development has been obtained from the analysis of simplified array 
models and from experimental data. The great speed and storage 
capacity of present day digital computers, however, have now made 
it possible to solve the planar phased array boundary value problem 
very accurately. 1 ' 2 

A general formulation of the planar phased array boundaiy value 
problem may be found. A vector two dimensional integral equation 



* The work reported in this paper was supported by the U. S. Army Materiel 
Command under contract DA-30-069-AMC-333(Y). 
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for the tangential aperture field (that is, the tangential field at the 
planar interface between the waveguides and free space) can then 
be derived. 

In its most general form the array elements are assumed to be 
arranged in a doubly periodic grid along two skewed (nonorthogonal) 
coordinates; and dielectric loading, covers, and plugs, as well as thin 
irises at the aperture plane, may be accounted for in the analysis. 
The possibility of multimode excitation of the array has also been 
included. The Ritz-Galerkin method is applied to obtain a solution 
for circular waveguide arrays. 

Numerical solutions for the reflection characteristics of dielectric- 
free planar arrays of circular waveguides hexagonally arranged in a 
conducting ground plane have been carried out. Experimental meas- 
urements have been made which compare favorably with the results. 
Forced surface waves 3-11 are found to occur at isolated points in the 
scan coordinates and can be related to certain vector and geometrical 
symmetries for an equilateral triangle grid array. These surface waves 
(or resonances) are often difficult to locate experimentally by the use 
of waveguide simulators 12-14 or small test arrays. The strong depend- 
ence of these forced resonances upon the mode of excitation is also 
demonstrated. The reflection characteristics as well as the polarization 
characteristics of the radiation pattern are illustrated for various 
combinations of linear and circular polarization excitation of the 
array. 

II. ANALYSIS 

An infinite planar array of waveguide elements, Fig. 1, is imbedded 
in a conducting ground plane at its interface (plane z = 0) with free 
space. The elements are arranged in a periodic grid along the skewed 
(nonorthogonal) coordinates s a and s 2 . The x and s 2 axes coincide 
while the s 2 axis makes an angle a with respect to the x (and s x ) 
axis. The element location is defined by two indices (p, q) correspond- 
ing to a physical location 

p po = ptei + qds 2 (1) 

where Si and s 2 are unit vectors along the Si and s 2 axes, while b and d 
represent the basic periods of the two dimensional grid. A basic 
periodic cell, 15 the parallelogram shown in Fig. 1, is thereby defined. 
If the array elements are excited uniformly in amplitude with a linear 
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Fig. 1 — Circular waveguide array geometry, 
phase taper such that the phase of the (p, q)th element is 

hi = P^ii + #«, , (2) 

then the resulting electromagnetic fields in the {p, q)th and (m, n)th 
periodic cells satisfy the following periodicity relationship 

E( QpQ ) = E( Pmn ) exp [i {(m - p)+ tl + (n - g)*J] (3) 

where E(p po ) may designate the electric or magnetic field at the (p, q)th 
periodic cell of the grid.* Therefore, except for a phase factor, the fields 
in all the cells are identical. 

In order to solve the boundary value problem, the exterior (free space) 
fields are expressed in terms of a complete set of Floquet type solutions 
of Maxwell's equations { *F* n exp ± iB mn z). These vector modes, which 
are functions of the steering phases \{/ tl and \f/ t , are derived in Appendix 
A. The interior (z ^ 0) fields are expressed in terms of the appropriate 
waveguide complete orthonormal set of vector modes {«!>/ exp ± «T,z} .f 
The boundary value problem is expressed in terms of an integral equa- 
tion which includes the necessary continuity conditions. This equation 
is formulated 1 by satisfying the continuity of the transverse (to z) 



*The parallelogram in Fig. 1 defines the (p = 0, q = 0) cell. The (p, g)th cell 
is translated by pb and qd along the Si and s 2 axes, respectively. 

t The waveguide modes are real functions and in general consist of both TE 
and TM modes with double subscripts. However, one can always systematically 
relabel these modes with a single subscript according to the increasing values of 
the eigenvalues. 
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electric and magnetic fields within a single periodic cell. As shown in 
Appendix B, the periodic cell consisting of the parallelogram CDEF of 
Fig. 2 can be replaced, without a loss of generality, by the parallelogram 
GHIJ (or any other periodic contiguous cell containing a complete 
single waveguide aperture). 




G D H 

Fig. 2 — Periodic cell in skewed arraj' geometry. 



The tangential electric and magnetic fields at the array interface 
(E and H at z = 0) can be expressed in terms of a Fourier series of the 
complete orthonormal set of waveguide modes {<*>, } for z ^ and by the 
set of Floquet type modes {W v mn ) for z ^ 0. Let the waveguides be 
excited by any linear combination of their propagating modes* with 
amplitudes A t (j = 1, ... , J for J propagating modes) and let the 
coefficients R { represent the amplitudes of the corresponding reflected 



*It is straightforward to include, if desired, any linear combination of both 
propagating and evanescent modes. 
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modes.* Let the coefficients D, represent the amplitudes of the reflected 
evanescent modes which are generated at the aperture. Then, in terms 
of the waveguide modes, the electric field at z = is given by 

f J °° 

2 (A; -f- R f )&j + X) A-*/ ( over tne waveguide aperture) 

e_ = j'- 1 p \ ' j mJ y (4) 

(over'therest'of'theperiodic cell). 
The corresponding magnetic field is 



-H_ = 



ef,(4,-b,)*,- e y^i*i 

(over the waveguide aperture) . . 
_0 (over the rest of the periodic cell). 

where the modal admittances {Y,} are real for propagating modes and 
pure imaginary for evanescent modes. These admittances are given by 10 

Y { = — for TE modes; F, = ^ for TM modes (6) 



toy. 



r, 



for an exp [—tut] time convention with the T, (the z propagation con- 
stant) being positive imaginary for evanscent modes. The tangential 
electric field at z = + , expressed in terms of the Floquet type modes, is 

2 

E + = 2£2 F mnp *I'* n (over the periodic cell) (7) 

p-l (m) (n) 

where the superscript p designates TE(p = 1) or TM(p = 2) modes. 
The magnetic field is correspondingly given by 

2 

— H+ = /—I 2—1 2—1 ^mnp-* mnp *mn \°/ 

p-l (m) (n) 

where the modal admittances Y' mnt> are given by 

y, m Bjsn . y, J*.. (Q) 

From the orthonorniality of the sets {*,} and {*F* n } it is clear that 

+ Rj for j ^ J ^ 

3>J 






* Actually {R } } can represent the reflection coefficients once {A,} are properly 
normalized. 
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where A is the waveguide aperture, and 

<E + , W p mn ) = ff E,.(rj* da = ff E + -(W v mn )* da = F mn „ (11) 

peri oil i o col I A 

where the symbol * denotes the complex conjugate. Notice that E+ 
vanishes on the conducting ground plane. To insure the continuity 
of the tangential fields across the aperture at z = 0, one requires 

E+ = E_ = E< over the aperture and the periodic cell (12) 

while 



H + = H_ as H, over the aperture only. 



(13) 



Using the various relations, (4) through (13), one obtains an integral 
equation for the tangential electric field E, at the array interface 



2 fM,y>,- £i>, ff^-E, 



da 



+ ZZE YLjKm ff 0*1)* -*i da. 

Vl (m) (n) JJ 



(14) 



Similarly, one can obtain 1 an integral equation for H, which is defined 
over the entire periodic cell. Under certain conditions 17 it is possible to 
interchange the order of summation and integration in (14) and thereby 
obtain the usual form of Fredholm integral equation of the first kind. 

A useful method of obtaining a solution to (14) is by the application 
of the Ritz-Galerkin method, 18 whereby the integral equation is reduced 
to a linear matrix equation. Substituting (4) in (14) for E, and taking the 
moments of (14) with respect to the set {*,•}, while using (10), leads to 
the following matrix equation 

r.^,1 \a, + b, 



YjAj 



= \K 



Aj +Rj 
D J+l 



(15) 



where \K\ is a square matrix with the {i, q) th element given by 
K - F, *.-„ +EZE TLJZJOZU . (16) 

p=l (m) (n) 
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In (16) 8,q is the kroenecker delta 



5., = 



1 if i = q 



if i 9± q 



and 



C 



,, - // 



W mn <Pi da 



(17) 



(18) 



is the coupling coefficient between the designated interior and exterior 

modes.* 

The practicality of an accurate numerical solution or approximate 
analytical solution often hinges upon obtaining C p mni in closed form. 
Recently the authors 10 have obtained closed form expressions of these 
coupling coefficients for circular as well as coaxial waveguides. 

A solution for the aperture field in terms of the coefficients of the 

waveguide modes is given by 



A, + R x 

Aj + Rj 
D J+l 



= 2\K 



A 1 Y J 

AjYj 





(19) 



The solution vector can be obtained by matrix inversion or by rapidly 
convergent iterative methods. 5,20 A similar procedure should be followed 
to obtain H, except that the aperture field and the moments may be 
taken with respect to the set { W mn \ . 

Once the aperture electric field is obtained, the input impedance 
and radiation properties of the array, as a function of scan, are easily 
obtained. The reflection coefficients are obtained directly from (19), 
as are the amplitudes of the evanescent modes in the waveguides. The 
radiation pattern of a single element in the array environment, includ- 
ing its polarization characteristics, can also be easily obtained. 21 - 22 

The addition of either a dielectric sheath or plug or both to the 

* Notice that other sets of functions { £j } can be used to reduce ( 14) by the method 
of moments. However, the integrability of 

/J u )*;•?' a nd Siu) x V^-U 
may prove difficult depending on ^i . The convergence properties as a function of 
the order of K will also be influenced by the choice of (&}. 
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array, Fig. 3, does not alter the functional form of the integral equa- 
tion (14). As was shown by Galindo and Wu, 3, 23 only the modal 
admittances in (14) have to be replaced. For the case when dielectric 
plugs are used, Y } is replaced by 

Y, -iy, tan yfr (m 

where y f and 7; are the modal admittance and propagation constants, 
respectively, of the <t>, mode in the dielectric and t x is the dielectric plug 
thickness. The relation between the reflection coefficients of the propa- 
gating modes in the dielectric-free region (e ) and the aperture field is 
given by 

y t J J E r <P, da - A,- exp(— iT/t^y, cos 7,^ + iT,- sin -y^] 

E> _ A 

exp(tT,0[y,- cos 7,*i - iYi sin 7^] ' 

the phase of Rj being referred to the aperture. Similarly, for the 
exterior dielectric sheath covering the array, the modal admittances 
are replaced by 

Y mnP - y«p yLp _ iYLp tan ^ , (22) 

where y' mnv and /3 mn are the modal admittance and propagation constants, 
respectively, of the W v mn mode in the sheath, and t[ is the sheath thick- 
ness. The coefficient of the mode *F* „ in the free space region above the 
sheath, F mnp , is related to the aperture field by the following relation: 

r-*-. 7 ( ,r g -ff y -T - ,, ff 1.-VS*. (23) 



y m ' np cos p mn ti — iY' mnv sin p mn t[ 



iodio 
■ 11 



the phase of F mnv being referred to the aperture plane. 

The integral equation formulation can be extended to the case when 
thin metallic irises are present at the waveguide aperture for matching 
purposes (see Fig. 4). The integral equation for the aperture electric 
field, (14), is still valid except that the integral has to be defined over the 
effective aperture with the result that the orthogonal relations of (10) 
and (11) cannot be used for the Ritz-Galerkin method of solution.* 
The modal coupling coefficients, (18), are still integrable in closed form 

* An integral equation for the magnetic field is not valid in this case because of 
the discontinuity of H< across the iris. 
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Fig. 3 — Dielectric sheath and plug geometry. 



over the effective aperture if the iris is circularly symmetric. However, 
the matrix elements, as given by (16), do change in form when an iris is 
present at the aperture. As a result of the <I>, not being orthogonal over 
the new effective aperture, the term F.-5,-, in (16) is replaced by an infi- 
nite sum. 

Multimode excitation of waveguide antenna fields has been used for 
primary pattern control. Such excitations may prove useful for the 
reduction of mutual coupling effects. 7 They also may be used to obtain 
circular or elliptical polarization from two linearly polarized modes and 
improve the polarization characteristics of the array. In the circular 
waveguide array, the horizontal TE n (*0 and vertical TE n (<E> 2 ) modes 
are degenerate (that is, they have the same z-directed propagation 
constant and impedance). In order to obtain a linear, elliptical or 
circular polarization excitation of the waveguide, one may redefine the 
first two modes as 4> liV and 4> 2 n' 



*P lN = 



A, 



(I A, | 2 + | A 2 | 2 ) 



T*i + 



(| A, | 2 + | A 2 | 2 ) 



-\* 



<l>, 



(24) 



EFFECTIVE 
APERTURE 



ru 



Fig. 4 — Aperture iris geometry. 
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*,v = - , , ,,'r, 1 , iW #» + Ii7i ,%',,', gff ^ • (25) 



A* | A 
'(l -A. | 2 + \A 2 l 2 )^ 1 ' A?(| Ax | 2 + | A 2 H 
This redefinition of the first two modes preserves the orthonormality 
and completeness of the set of waveguide modes while allowing the 
flexibility in adjusting the desired array excitation. The reflection 
coefficients which correspond to these redefined modes are R^ and 
R 2 n- The polarization characteristics of the radiation pattern may 
be determined from the and <f> components of the radiated field, 
Eg and E+ , respectively. After proper normalization of Ee and E^ one 
may obtain the corresponding transmission coefficients 19 



T s = 



■» (1(12 \ -•- 



l2\* 



YJ (1 - Tl - Tl) 

(aperture) 



fj E r (W 2 00 )*da (26) 



T ' = ~{ X t) i If E «"«o)* da (27) 

(aperture) 

where T„ and T v define the scan angle directional cosines. When the 
first two waveguide modes are the only propagating ones and while 
only a single lobe propagates in the free space, equations (24) through 
(27) are related by the conservation of energy relation: 

\Rxn\*+ \r 2N \ 2 + 1 r# r + 1 av r - !■ (28) 

For more than one lobe in free space or additional propagating modes 
in the waveguides, (28) has to be accordingly modified. 

III. NUMERICAL AND EXPERIMENTAL RESULTS 

In order to obtain a numerical solution for the aperture field, the 
infinite dimensional matrix of equation (15) must be truncated and 
cast in a finite dimensional form. In other words, the electromagnetic 
fields will be approximated by a finite Fourier series of the waveguide 
and free space modes, and consequently the solution of the problem 
as given by (19) , is finite dimensional as well. 

In numerical solutions of problems of this type, various ways of 
ascertaining the validity and accuracy of the solution are desirable. 
One obvious way is to increase the number of waveguide and free 
space modes and check the convergence of the solution as a function 
of the number of modes used in truncating (15). However, for the 
type of kernel involved in this problem, monotonicity of the conver- 
gence is not assured. Nevertheless, convergence is an important check 
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since this numerical solution is variational or stationary for the im- 
pedance. 12, 24 Iterative methods and error estimates 5 - 20 may also be 
used for checking the convergence of the solution. Special symmetries 
of the reflection coefficients versus the scan angles which are dictated 
by the array geometry and mode of excitation can serve as a semi- 
independent check. An important independent check in which the 
reflection and transmission coefficients can be measured at special 
scan angles with the aid of a waveguide simulator is used as well. 
The numerical results will in general be presented as a function 
of scan angle. For convenience, however, the differential phase shifts 
between elements will be used as the independent variables. These 
quantities are fa in the ^-direction and $ v in the ^/-direction. Further- 
more, since we limit ourselves to radial planes of scan, we introduce 
the quantity if/,-- These quantities are related to the directional cosines 
as follows: 

*- 2 r r - s ^ = M f^ T '- *-«s + rf>'- ( 29 ) 

The amount of computation can be reduced when one recognizes that 
the following symmetry in the aperture field as a function of scan 
exists : 

E,0fc,*,) = E,(-*., -*„). (30) 

Convergence tests as a function of the number of waveguide and 
free space modes indicated that 18 waveguide modes and 338 free 
space modes yield several percent (usually less than 2 percent) 
accuracy in the magnitude of the reflection coefficients, Rj, except 
near sharp changes of Rj where the position of the sharp changes is 
accurate to several degrees in ij/ r . 

The energy conservation check, equation (28), is a necessary check 
but not a sufficient one in this problem as well as in other interior 
type boundary value problems. 17, 22, 33 

One of various special symmetry checks is depicted in Fig. 5. A square 
grid array in the (x, y) coordinates is excited by the vertical TE U (4> 2 ) 
mode. In this coordinate system the array parameters are represented 
in the following way 

a = waveguide radius; b = d; a = 90°| ,„.,. 

®IN = *2 ; ®2N = — *l J 

The parameters of the same array, when viewed in the (x', y') coordi- 
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y 



d=b=d 




Fig. 5 — Symmetries in square grid arrays. 

nate system, can be alternatively represented as 



V = (2)*d' = (2)»d;"j 



p*' = 45° J 



1 



-*i +*2}; 



l 

w 



{&{ +<!»£: 



(32) 



where the reference to the (a/, y') system is denoted by primes. The 
same results should be found for this array at any scan direction 
regardless of the representation. This offers, to a degree, a check on 
roundoff error. Numerically, the reflection coefficients differed only 
by a fraction of a percent. 

An additional symmetry check is given in Fig. 6, where the magnitudes 
of the reflection coefficients are plotted versus \f/ r for a 45° plane of scan. 
Ri N and R 2tf correspond to the reflection coefficients of ® 1N and <S> 2 n as 
defined by (31). At fa = 240° (shown by a vertical arrow) the main beam 
is grazing, while for fa > 240° no beam exists in real space and the total 
incident power is reflected and divided between the two propagating 
modes, «*»!# and «*> 2 ;v . Of special interest is the point fa = 180° 
X (2)* pa 255°. At this point fa = ft, = ^,, = fa, = 180° and the array 
excitation is as indicated in the inset of Fig. 6. If the array is to be 
simulated at this scan angle, the appropriate waveguide simulator would 
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144 180 216 252 28{ 

f r IN DEGREES |fl = 90° 



Fig. 6 — Square grid array: \Rur\ and |/?iv| vs t/- r in the 45° plane of scan. 
a = 0.48, b = 1, d = 1, \ = 1.5 and a = 90°. 



consist of the square waveguide (solid lines) shown in the inset. It is 
clear that the horizontal waveguide mode (*P 2 n) cannot be excited. The 
numerical results indeed show that | R 2 n | = at this scan angle. 

Figure 7 shows a close agreement between experimental and numeri- 
cal results for a rectangular grid array with vertical polarization ex- 
citation scanned in the H-plane. The scan angle which corresponds to 
steering phases \f/ x = \f/ tl = 180° and yp y = \f/ tl = 0, can be simulated 
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Fig. 7 — Rectangular grid array: computed and experimental results vs fre- 
quency at //-plane symmetry point, a = 3.57 cm., b = 16.51 cm., d = 8561 cm., 
and a = 90°. 



1916 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1968 






Fig. 8 — Rectangular grid array: simulator for .ff -plane symmetry point. 

by the rectangular waveguide (solid line) in the inset. The experi- 
mental results were obtained from measurements of an abrupt junc- 
tion between a circular waveguide and an L-band rectangular wave- 
guide as shown in Fig. 8. 



IV. EQUILATERAL TRIANGLE GRID ARRAYS 

Let us consider the reflection and radiation characteristics of cir- 
cular waveguide arrays arranged in an equilateral triangle grid, and 
the strong dependence of the array properties upon the mode of 
excitation. Grating lobe incipience or a beam at grazing is designated 
by a vertical arrow in the illustrations, 
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Figure 9 shows the reflection coefficient of an array in the E-plane 
scan, with vertical polarization excitation. In this plane the hori- 
zontal mode is not excited because of symmetry so that R 2 n = 0. As 
can be seen, the slope of both the magnitude and phase of R iN is dis- 
continuous and singular at a grating lobe incipience, which parallels 
previous observations in rectangular waveguide arrays. 26 - 27 This is 
related to the asymptotic decay of the coupling coefficients. A forced 
surface wave resonance can be seen around grating lobe incipience 
where | Rw | = 1.0. Notice that this forced surface wave resonance 
is extremely sharp and consequently may not be observed in small 
finite arrays. 

The corresponding transmission coefficient Tg is shown in Fig. 10. 
The plot of the transmission coefficient is actually the radiation pat- 
tern of a single element in the array environment and it exhibits the 
null which corresponds to a total reflection. Notice that the phase 
of the transmission coefficient will exhibit a 180° discontinuity when 
the magnitude has a zero. The magnitude of the reflection coefficients 
of the same array, in a 60° plane of scan, are shown in Fig. 11. Again, 




Fig. 9 — E-p\ane scan of Rw vs <^ r (a = 0.48, 
b = 1, d = 1, X = 1.4 and a — 60°). 
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fp IN DEGREES Iff = 53.9" 



Fig. 10 — ff-plane scan of T vs f r (o = 0.48, b = 1, d = 1, X = 1.4, and « = 60°). 




144 180 216 

iK, IN DEGREES 



252 288 



360 



Fig. 11 — \Rin\ and \Rts\ vs i/- r in the 60° plane of scan (a = 0.48, b = 1, d = 1, 
\ = 1.4, and a = 60°). 
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the singular slope at grating lobe incipience can be observed. As can 
be seen, the distribution of the reflected power between the two modes 
is a function of the scan angle. 

An interesting phenomenon can be observed in the corresponding 
transmission characteristics which are shown in Fig. 12. The <£ trans- 



i - 



^ 



c u> </) a 





Fig. 12 — \T \ and |^| vs ^ r in the 60° plane of scan (a = 0.48, b = 1, d = 1, 
\ = 1.4, and « = 60°). 

mission coefficient, T^, vanishes prior to grating lobe incipience (the 
positions of the vertical arrows). The vanishing of T$ at this scan 
angle can be directly related, when coupled with the vector symmetries 
in the array excitation and geometry, to a forced surface wave reson- 
ance. If the array excitation consists of the sum of the two modes 
(equal in phase and magnitude) indicated by the solid arrows in 
Fig. 13, then in the 60° scan plane To = by symmetry considerations 
and T, t , vanishes as shown in Fig. 12. Since the vectorial sum of the 
two solid arrows in Fig. 13 is the dashed arrow, zero transmission or 
a forced surface wave resonance will occur in the /f-plane of this 
polarization. Figs. 14 and 15 indeed show this effect in the i?-plane 
scan where | R iy | and | T^, | attain unity and zero respectively prior 
to grating lobe incipience.* 

Since the forced surface waves are related to the vector symmetries 
just mentioned, one may anticipate that the scan points at which they 
occur are isolated. Figure 16 indeed demonstrates that the forced reso- 
nances in the E and H planes occur at isolated points. The scan around 



* The difference in the values of \p r at which these phenomena occur is inherent 
in the definition of \f/ r , equation (20) . 
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Fig. 13 — Vector symmetry relationships between transmission zeros and forced 
surface waves. 




Fig. 14 — //-plane scan of Rik vs $ r (a = 0.48, b = 1, d = 1, X = 1.4, and a = 60°) 
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Fig. 15 — #-plane scan of T« vs \f/ r (a = 0.48, b = l,d = l,\ = 1.4, and a = 60°). 

the grating lobe circles in Fig. 17 shows that the peak of the total reflec- 
tion \ R T \ = (|^iat| 2 + I R2N 1 2 ) varies from unity in the H plane 
(8 = 0), gradually decreases, and then increases again and reaches unity 
in the Z?-plane (8 = 90°). The anomalous behavior near grating lobe 
incipience is eliminated when the polarization of the excitation is changed 
to horizontal, as shown in Fig. 17, indicating thus the strong dependence 
of the forced surface wave resonances upon the mode of excitation. Over 
the frequency band given by 1.3 ^ X ^ 1.5, qualitatively similar 
behavior of the radiation and reflection characteristics of the array was 

observed. 

Figure 18 shows the reflection characteristics of the array under 
circular polarization excitation. In this case 



O lV = — 



1 

w 



<!»! + »*2j 



Oov = — 



(2)' 



4»! — Z4» 2 



(33) 



The incident mode is 4> IA r. 

Again the singular slope of these curves at grating lobe incipience 

can be seen. The division of the reflected power between the two 

modes as a function of scan may be observed as well. Fig. 19 shows 

the polarization characteristics of the radiation pattern of an array 

element. The axial ratio, denoted as A.R., is the ratio of the minor 

to major axis of the polarization ellipse while the tilt angle of the 

major axis, t, is taken with respect to the <£ axis. As indicated by 
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Fig. 16 — |i2r| vs \j/ r in the vicinity of grating lobe incipience. 
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Fig. 17 — E- and ff-plane scans of |it!iAr| vs ^ r for horizontal polarization (a = 
0.48, b = 1, rf = 1, X = 1.4, and a = 60°). 



CIRCULAR WAVEGUIDE ARRAYS 



1923 




144. 180 216 

i/i IN DEGREES 



Fig. 18— \Rxa\ and \Rw\ vs \p r in the 0° plane of scan. Circularly polarized ex- 
citation (a = 0.48, b = 1, d = 1, X = 1.4, and a = 60°). 

the plot, the element (or the array) far field pattern is circularly 
polarized around broadside, A.R. = 1.0. It deteriorates to linear 
polarization, A.R. = 0, at two points prior to grating lobe incipience. 
The linear polarization at $ r = 203.5° with t = 90° results from the 
//-plane forced surface wave resonance of Figs. 14 and 15 where 
r vanishes. The null of the axial ratio at ^,. = 207° (with t = 105°) 
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Fig. 19 — Radiated axial ratio (A.R.) and tilt angle (t) vs fr in the 0° planr 
of scan. Circularly polarized excitation (a = 0.48, b = 1, d = 1, X = 1.4, and 
a = 60°). 
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is caused by the difference between the phases of T g and T^ causing 
the 4> and 8 components of the far field to be in phase. 

The reflection and polarization characteristics in the 30° plane of 
scan are shown in Figs. 20 and 21, respectively. In this case the 
single null of the axial ratio at grating lobe incipience results from 
the Z?-plane forced surface wave resonance of Figs. 9 and 10. For 
planes of scan between 0° and 30° the results corresponding to one 
plane change gradually to those of the other plane. Around grating 
lobe incipience the axial ratio drops appreciably (to around 0.1) but 
does not reach zero. From the circular symmetry of excitation and 
six-fold symmetry of the array geometry, a 30° sector of scan com- 
pletely specifies the array reflection and radiation characteristics. 

V. CONCLUSIONS 

A general formulation of the planar phased array boundary value 
problem has been given in terms of a vector two dimensional integral 
equation. The solution of this equation by the Ritz-Galerkin method 
closely agreed with experimental results. 

Equilateral triangle phased arrays of circular waveguides were 
numerically analyzed. It was found that forced aperture resonances 
or forced surface waves, manifested by total reflection and no radia- 
tion, do exist for these arrays even in the absence of dielectric ma- 
terials. These effects were observed over a 15 per cent frequency 
band. The forced aperture resonances occurred prior or close to grat- 
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Fig. 20 — |i?i-v| and |fl».v| vs \p r in the 30° plane of scan. Circularly polarized 
excitation (a = 0.48, b = 1, d = LA = 1.4, and a = 60°). 



CIRCULAR WAVEGUIDE ARRAYS 



1925 




Fig. 21 — Radiated axial ratio (A.R.) and tilt angle (t) vs <p T in the 30° plane 
of scan. Circularly polarized excitation (a = 0.48, b = 1, d = 1, X = 1.4, and 
a = 60°). 

ing lobe incipience in the E and H plane of scan for vertical polariza- 
tion excitation. These resonances were found to occur at isolated 
points as a function of the scan variables and are strongly influenced 
by the mode of excitation. The resonances vanish when the polariza- 
tion of excitation changes from vertical to horizontal. 

The polarization characteristics of the radiation pattern (or al- 
ternatively the radiation pattern of a single excited element in the 
array environment) is shown at selected planes of scan for circular 
polarization excitation. The degradation of the axial ratio resulting 
from the forced surface waves was shown. Total reflection or no 
transmission owing to forced aperture resonances were not observed 
for circular polarization excitation in the cases presented. 

The analysis of coaxial waveguide arrays, as well as the incorpora- 
tion of thin, circularly symmetric irises in the aperture of the wave- 
guide clement, can be carried out along lines similar to those dis- 
cussed here. 

The effects of dielectric loading of the array as well as dielectric 
covers and plugs have also been studied. 34 
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APPENDIX A 

Floquet Type Wave Functions in Skewed Coordinates 

Consider the periodic array, Figs. 1 and 2, excited with incremental 
phase shifts (steering phases) between adjacent cells along the s x and 
s-2 coordinates. A complete set of solutions to the scalar wave equation, 
each of which varies periodically, according to Floquet's theorem, 
along the «i and s 2 coordinates is 



S mn = exp (iB mn z) exp i 



^j, 2irm\ 



.[>,, 2™"] 



s t exp i\—r j- s 2 (34) 



with the integers m,n= — =0 , . . . — 1, 0, 1, 2, 3, . . . , + °° • Equation 
(34) describes a wave traveling (or decaying) in the 2 direction with 
propagation constant B mn (exp — iwt time convention). The steering 
phases \f/ tl and ^,, are directly related to the beam pointing direction, f, 
of a radiated plane wave with a vector propagation constant k = k Q f, 
so that (34) can be rewritten as 

S mn = exp (iB mn z) exp i\ ko-Sj - -\- [s, exp i k -s 2 - —3- \s 2 . (35) 



b I 



d J 



The free space propagation vector k can be expressed in the cartesian 
coordinate system as 

k = k [TJ 4- T v y + TJ] (36) 

where T z , T u and T t are the directional cosines of k with respect to 
that system and x, y, and z are the unit vectors. The quantities 

ko-Si = ~ and k -s 2 = ^ (37) 

are the projections of k on the reciprocal grid (lattice) coordinates t t 
and t 2 , respectively. 15,28-30 The unit vectors in the ti and t 2 directions 
form a biorthogonal set with s, and s 2 (Fig. 22). To express (35) in 
cartesian coordinates it can easily be shown that 

Sj = x — y cot a, s 2 = y/sin a, (38) 
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Fig. 22 — Grating lobe diagram in reciprocal lattice coordinates. 

so that the substitution of (38) into (35) yields 
S mn = exp (iB mn z) exp i\ k T, - —j^r 



.exp{fc o r i/ -(^^- 



2irm 
b tan i 



y. (39) 



The propagation constants of the {m, n)th Floquet mode along the 
x and y directions, k x and k„ respectively, are 



K z = K • X — KqI z 

k u ^k-y = A- T„ - 



2irm 



b ' 
2tm 



(40) 



2-irm 



d sin a b tan a 



Since S mn is a solution to the scalar wave equation, it can be shown 
that 
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B mn = (*" - kl - k 2 rf 



= \K - 



k T,-*f\ 



■I-[ 



k T u — 



2ir?l 

d sin < 



J,wm 
b tan < 



(41) 
1 



where the positive imaginary root holds for (k 2 x + K) > h% (time 
convention exp — iiot). 

Each mode, S mn , for which B mn is real corresponds to a radiated plane 
wave of the phased array. The plane wave with the indices m = and 
n = is identified with the main beam while m or n ^ corresponds to 
a radiating grating lobe. As a function of T x and T„ (or \f/,, and ^,J, a 
given B mn may become pure imaginary as it goes through a zero, as in 
equation (41). In such cases, the related Floquet mode, S mn , becomes 
evanescent or nonradiating. By plotting the curves obtained by setting 
the B mn = as a function of T x and T v , one obtains a convenient dia- 
gram which illustrates these effects. Setting B mn = yields 



[,. _ ?f] + [,. _ ( 



n\ 



mX 



d sin a b tan a 



= 1, 



(42) 



where Ao = 2ir/k . As a function of T„ and T v , (42) represents a 
family of circles with unit radius displaced from the origin. This 
diagram of displaced circles constitutes the well-known grating lobe 
diagram, Fig. 22. 31 - 32 

Notice that the steering phases, \f/ tl and \fj t z, are related to T x and 
T v through equations (36) and (37) . The parallelogram C'D'E'F' of 
Fig. 22 corresponds to the range of steering phases 



— r ^ ^ tx ^ ■*; 



- » 25 *«. £ 



(43) 



and is a periodic cell along the U and t 2 coordinates. 

As mentioned in Section II, it is possible to define a complete ortho- 
normal set 16 of vector modes {^*^ n } over the parallelogram CDEF, 
Fig. 2. The tangential electromagnetic field at the plane z = + can be 
expressed by a Fourier series of this set of modes which consists of both 
TE and TM modes (transverse to z). These modes, { W T m E n } and { »F™ } , 
are given by 



and 



^ = 



* win 



exp i(xk z + yk u ) Jfc 
(bd sin a) h 



K 



kr X ~k r y 



exp i(xk x -f- yk u ) 
(bdsino:)* 






(44) 



(45) 
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with K = (k\ + K)K The quantities k x and k u are functions of (m, n) as 
given by (40). The orthonormality of this set of vector modes is defined 
by the following scalar products: 



■;•)= // vj.i 



(vS , <,■> = J J v2 • « E )* <fc rfy - *...p. , (46) 

where 



pa ml le locram 
CDEF 



and 



APPENDIX B 



[l for m = p and n = q, ,^ 

[0 otherwise 
OrS , *F™> = (48) 



On the Invariance of the Scalar Product with 
the Shape of the Periodic Cell 

The orthonormality and completeness of the set {W v mn } of Floquet 
modes, equations (44) and (45), need not be denned over a specific 
periodic cell such as the parallelogram CDEF of Fig. 2. This fact is 
especially significant when a periodic cell intercepts parts of more than 
one circular (or other type of element) aperture. It will be shown that 
the orthonormality of { W p mn ] can be preserved over a properly deformed 
periodic cell which contains only a single waveguide aperture. 

Using (37), equation (44) can be rewritten as 

Wll = F(w?, n) exp i\ ^— ^ s, + ^^ s 2 I . (50) 

The scalar product betewen two TE modes is 

\**mn > * VQ I "mn.pq 



= // ^(<, E )*sin<*d Sl ds 2 

(CDEF) 

= F(m, n) -F(pg)* sin a 



(51) 
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// exp -i[ Mm b - P) ., + ^^ *,] A, A. ■ 



The integral over the parallelogram CDEF in (51) can be divided 
into three (or more) integrals over the triangles CGK and LEI and 
the polygon GDLIFK (Fig. 2) : 



CDEF CGK LEI GDLIFK 



(52) 



Because of the array periodicity, the triangles DHL and LEI are 
displaced by b with respect to the triangles CGK and KFJ, respectively, 
along the s, direction. Thus, for example, if the s, coordinate of the points 
within the triangle DHL is s[ given by 



s[ = s, + b with ds[ = ds x 



(53) 



then 



// exp -(M^ s ; + ?^p) *] rfs ; * 



= exp —i2ir(m — p) 1 1 exp — i\ r — ^ s, H ■= — — s 2 ds t 



= // exp -i 



dsi 



'^g) si + 2x(n- g ) g j^ 4i 



(54) 



Similarly 



dfia 



" // GXP " 



f 2 T (m - p) , 27r(n - q) 
1 6 Sl+ d 



uS, ds 2 . 



(55) 



Thus 



CDJ?/' 

" K GXP ~ { 



Mm-v) Si + MLjL& 82 ] d8i 



ds 2 (56) 
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and the orthonormality and completeness of the set { *F*»} is preserved 
over the new periodic cell GHIJ. This is the two dimensional analog of 
the single dimensional Fourier series whereby the functions and coeffi- 
cients are independent of the initial value of the period. In fact, the two 
dimensional periodic cell can be deformed to any shape which contains 
a single waveguide provided that the area of the cell stays the same and 
the parts of the cell which cause the deformation are translated by b or d 
along the s, and s 2 coordinates, respectively. 
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